We are IntechOpen, 
the world’s leading publisher of 


Open Access books 
Built by scientists, for scientists 


9,200 128,000 150M 


ailable International authors and editors Downloads 


Our author among the 


104 TOP 1% 12.2% 


Countries delivered to most cited s Contributors from top 500 universities 





Selection of our books indexed in the Book Citation Index 
in Web of Science™ Core Collection (BKCI) 


Interested in publishing with us? 
Contact book.department@intechopen.com 


Numbers displayed above are based on latest data collected. 
F information visit www.intechopen.com 


y 


5 


Several Kinds of Modified SPIHT Codec 


Wenchao Zhang 
Institute of Electronics, Chinese academy of science 
People’s Republic of China 


1. Introduction 


SPIHT (Set Partitioning In hierachical trees), being an efficient coding method for wavelet 
coefficients, has acquired more and more widely application, especially in image/video 
compression fields. But, conventional SPHIT have some obvious limitions. For example, 
when for the color image compression, polarmetric SAR image compression, or multi- 
spectrum image compression and other multi-channel image compression, there are only 
very limited image planes(R,G,B for color image, HH, HV,VVfor polarmetric SAR images) 
but there exist large amount of information redundancy among the image planes. So, 
considering the support length of discrete wavelet transform, we can’t use the set 8-partition 
methods such as 3D-SPIHT which has been used for video compression. Another example, 
when the input image is unsymmetrical such as 16x512 image block even only one line 
image content, which is very common in hardware design, because it means the larger final 
chip die size for the many line buffers. But, the traditional SPIHT codec can acquire the best 
compression performances only for the image is symmetrical in horizontal and vertical 
dimension. To the above specific applications, I will discuss several kinds of modified 
SPIHT in this chapter, most of them are the author's newly research result. 

In the following section, We will give a simple description about the traditional SPIHT 
codec, then, we will take the polarimetric SAR intensity image compression as example and 
give some specific compression method for multi-channel image compression. Certainly, 
before encoding for the wavelet coefficients, we need a 3D matrix transform to remove the 
information redundancy among the image plane and in each image plane. For the 
unsymmetrical image compression used in hardware design, an unsymmetrical SPIHT 
codec is detailed addressed, at the time, its specific case for only 1 line image compression, 
1D SPIHT codec is also given. 


2. Traditional SPIHT codec[1] 


SPIHT can be used to compress the traditional image, but it need 2D discrete wavelet 
transform (DWT) before encode the wavelet coefficients. The correlation of DWT coefficients 
not only exist in each subband inside but also among different subbands. Conventional 
SPIHT codec is to encode coefficients by SOT (Spatial Orientation Tree), which utilize the 
correlation of the same subband and different subbands. Conventional SPIHT codec, having 
many advantages such as simpleness, embedded code stream when compared with other 
encoders, is a very efficient compression method. 
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According the SOT structure, the set partitioning can be defined as: 


T(i,j) = c(i,j) + D(i, j) 
D(i,j) = O(i, j) + L(i, j) (1) 
L(i,j)=>D(k,l) (k,1)€O(i j) 


Here, T(i,j) is spatial orientation tree and c(i,j) is a root node of the tree; D(i, j), O(i,j) 
and L(i,j) are node c(i,j)’s descendant node set, direct descendant node set and indirect 
descendant node set. Direct node set O(i,j) can be further partitioned into 4 nodes just as 
the following: 


O(i, j) = (c(2i, 25), ci, 2j * 1,c Qi * 1, 2j), c(2i+1, 2j 1)] (2) 


Conventional SPIHT encoding process can be divided into sorting pass and refinement pass. 
During the encoding process, 3 lists are used to record the corresponding encoding 
information, which include List of insignificant set (LIS), list of insignificant pixel (LIP) and 
list of significant pixel (LSP). The basic operation is significance test and the significance test 
function is just as the following: 





1; max C; j } 22" 
eu (T) = (i,j)eT (3) 
0, otherwise, 


The encoding process can be simply described as: 

Output the initialized threshold n = | log; (max|c(i, j) | ; set LSP as NULL; add the root nodes 
having and no having descendant nodes to LIS and LIP respectively. 

For a spatial orientation treeT(i,j), its nodes are partitioned into nodes c(i,j) and 
descendant nodes D(i,j) according to formula 1. If node c(i,j) is significant, move it from 
LIP into LSP. If D(i,j) is significant, then partition D(i,j) into O(i,j) and L(i, j). The 
significant nodes of O(i,j) will also be moved to LSP according to the significance test. If 
L(i,j) is significant, then L(i,j) will also be partitioned 4 set and test the significance of 
every set. This process will be continued on. 

Sorting pass is for every element c(i,j) of LSP; output the n” efficient bit of c(i,j) for the 
current threshold n. 

Decrease the n and continue the sorting pass and refinement pass until meet the required bit 
number or the end of encoding. 

The decoding process is same to the encoding process, which, compared with EZW, make it 
un-necessary to transmit position information because the position information is 
embedded in the encoding process. 





3. Multi-channel image compression 


SPIHT, utilizing the information redundancy of 2D wavelet coefficients and adopting set 4 
division, can be used to compress 2D image data. Naturally, it can be enlongated to 3D 
video compression, that is to say, if we adopt 3D DWT to remove the information 
redundancy of video data and choose set 8-partition instead of set 4-partition to encode 3D 
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wavelet transform coefficients, this method is called 3D-SPIHT. But at some cases, we only 
need to compress multi-channel image, such as color images (R,G, B), polarimetric SAR 
image (HH,HV,VV), multi-spectrum image, ect. Because the third dimension usually have 
only limited image planes and can’t be processed by supported discrete wavelet transform, 
however there exist much information redundancy among each image channel. In this case, 
consider it’s simplity, 2D DWT for each image plane and 1D DCT can be used to remove the 
information redundancy. In the following section 3.1-3, taking polarimetric SAR image as 
example, 3D matrix transform and the related compression method including bit allocation 
based encoding method and 3D SPIHT Embedded method will be addressed. 


3.1 3D matrix transform for multi-channel image 

At first, 3D-matrix [2-3] is adopted to represent the multi-polarimetric SAR intensity images. 
The 1st, 27d and 3rd planes represent the HH, HV and VV polarimetric SAR intensity image 
respectively. Assuming the image dimensions are: MxN , the multi-polarimetric SAR 
images can be written as a MxNx3 matrix f(x,y,z) (x=0,1,2,---M-1,; 
y=0,1,2---N-1;z=0,1,2). In order to remove the redundancy of the matrix, DWT, KLT, 
DCT and other linear transforms can be used. 

Because there are only 3 components in polarimetric channel and KLT is dependent on the 
statistic properties of data, DCT other than DWT or KLT or other linear transforms is chosen 
to remove the redundancy of polarimetric channels. According to DCT transform definition, 
the DCT transform can be represented as: 


1 2 
F(x,y,0)=—= >" f(x/y,2) 

a (4) 
F(x,y,Z) = EE fey acos = (Z=1,2) 
z=0 


The 3D-matrix can be composed in other orders, such as HH, VV, HV acted as the 1st, 2nd 
and 3rd planes respectively. The components of like-polarimetric(HH and VV) have strong 
correlation, while the components of cross-polarimetric (HH/VV and HV) have weak 
correlation. Both the DCT theory and experiment results prove that the DCT coefficients of 
the matrix composed of HH, HV and VV are the most concentrated and that the decoded 
images have the least loss at the same coding rate 

After 1D-DCT transform, the data power of the 3D-matrix is concentrated into the 1st plane 
and the redundancy among three data planes decreases greatly. In order to remove the 
redundancy in every data plane of the 3D matrix, 2D-DWT is chosen. According to the 
definition of 2D-DWT transform, the image data can be decomposed into horizontal, 
vertical, diagonal and low frequency components after horizontal and vertical filtering. The 
low frequency component can be decomposed further. After many level discrete wavelet 
transform, the data power is concentrated to the low frequency components. After 1D-DCT 
transform and 2D-DWT transform, the power of the whole 3D-matrix is concentrated onto 
the top left corner. 3 level wavelet transform is adopted, so each data plane is decomposed 
into 10 subbands including LL, HLs, LH; HH;, HL», LH», HH2, HLı, LH;, HHı. Of all the 
subbands, the LL; subband of the 1st plane has the highest power. 

The 3D-matrix transform of Multi-polarimetric SAR intensity images (1D-DCT among 
polarimetric channel and 2D-DWT in each polarimetric plane) is illustrated in Figure.1. The 
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1D-DCT and 2D-DWT are linear transform, so, the operation 1DCT, and 2DWT,,y can be 
inverted in sequence, that is: 


IDCT, QDWT, „(f (x,y,z))) = 2DWT, ,(ADCT;(f (x, y,z))) (5) 


m-3 m2 m=1 O 


hai 


The 1st plane HH 


The 2nd plane HV 


The 3rd plane VV 


Fig. 1. Illustration of 3D-matrix transform of multi-polarimetric SAR intensity images. 


3.2 Bit allocation based on differential entropy encoding 

The three mixed coefficient planes have different energy, which means different importance. 
Different bit allocation scheme will greatly affect the quality of decoding images even at the 
same mean quantization bit rate. Our aim is to find the optimal quantization bits allocation 
to make the decoding images have the least distortion, which is a constrained optimal 
problem. The rate distortion function R(D) is unknown in most cases, so its solutions have 
become a hot research issue. Two variant iterate search [4], Dynamic programming [5] and 
R/D curve modeling [6] are proposed, but these methods yet haven’t been adopted for its 
enormous computation complexity or the local optimal not the whole optimal bit allocation 
searched. In this subsection, a new bit allocation method based on differential entropy is 
proposed. Compared with the existing bit allocation methods, it is easy to find a better bit 
allocation in whole with the proposed method. At the same time, the computation is not 
very complex. 

Let R,, R;, R} be the optimal bit rate for the three mixed coefficient plane respectively and 
Rr be the mean bit rate; Let D(R) be distortion rate function. So, the optimal bit allocation 
problem can be represented as: 


nn, (D,(R,) + Dy (R5) + D3(R3)) 
1752/53 (6) 
subject to R, +R, +R, =3Rr 
Adopt the Lagrangian multiplier, 
F(R,, R5, R3) = Dy (Ry) + D(R5) + D4(R5) + AR — Ry = Ry = R3) (7) 


www.intechopen.com 


Several Kinds of Modified SPIHT Codec 


71 





OF | 
OR, 
OF 
OR, 
OF 
ORs 


From 





we can acquire 





DR). , 
eR, 

DR) |, 
OR, 

€D3(Rs) _, 
ôR; 





(8) 


According to the inequation of rate distortion function for non-Gauss continuous 


information source, we have: 


h(U)- slog 2zeD < R(D) < slog 


(9) 


where h(U) is the differential entropy of information source and o? is its mean square 


error. 


At high bit rate, the lower bound of the inequation approaches the real rate distortion 
function for most probability distributed information source. So, we can let the lower bound 


equal to the real rate distortion function and then acquire: 


R(D) = h(U) - slog 2reD 


Further, we can acquire 
and 


For every mixed coefficient plane: 


R, =h(U,) - Zlog(-e2) 





R, = h(U,) - log(-zeA) 
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According to the constrain condition R, + R - R5 -3R,, we can acquire the final bit 
allocation for every coefficient plane: 











R, = + ZH Hl) Hl) 
Rs — Re d 2h(U;) CA —h(U3) (14) 
Ry Ry + Aad) 





From formula (14), we can acquire the optimal bit allocation for the three mixed coefficient 
planes, then the conventional SPIHT can be adopted to encode the coefficients for every 
plane according to the corresponding allocated bit rate. 


3.3 3D-SPIHT embedded encoding 

The formula (14) provides a method to compute the bit rate allocation, but the bit allocation 
accuracy strongly depends on the degree of the lower bound of formula (9) approaches the 
real rate distortion function. So, in most cases, we only can acquire an approximating 
optimal bit allocation. In this subsection, a new method named 3D-SPIHT embedded 
algorithm is proposed, which extends the conventional SPIHT algorithm to 3D case. It 
avoids bit allocation by adopting 3D-SPIHT to encode the three mixed coefficient plane 
entirely, so the optimal bit rate allocation can be ensured. 

During the 3D-SPIHT encoding process, the following sets and lists will be used: 

Poss) s Diytane (0, J) s Ou (1, J) and Laune (5 J) represent the iplane" 
root node, descendant's nodes of node(ij) offspring's nodes of node(ij) and indirect 
offspring’s nodes of node(ij) respectively. LISiptane, LIPipane and LSPipiane represent the 


coefficient plane's 


iplane" plane's list of insignificant pixel sets, list of insignificant pixels and list of significant 
pixels respectively, where iplane=1,2,3. 

Just as the conventional SPIHT algorithm, the set splitting procedure for every coefficient 
plane can be defined as following: 


Tiptane (tJ) = H ;piane (tJ) + Diptane (Er j) 
D plane (i, j) = Oiane (i, j) + Lone, j) (15) 
Ligne, J) = I, Diptane(k/)) (k,l) e Ovtane(t/ j) 


The 3D-SPIHT process can be divided into sorting pass and refinement pass, whose basic 
Operations is also significance test of set just as the conventional SPIHT algorithm. That is, 


1, maxi|C,, (Lj)(22 
iplane 
S, (Tui) = (1,7) Tiptane (16) 


0, otherwise 


where C;;,,(i,j) is the wavelet coefficient of coordinate (i,j) in the iplane" mixed 


coefficient plane. 
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The followings are the coding process: 
1. Initialization 


n_max, = |1082 (maxi, (CiD) 


n_max, = |1083 (maxi, (el) 


n_max, = |108 (max¢, sell) 


n =n_max, 


Output (17) 








(For the HH HV VV combination, nr. max, » _max3 >n_max, can always be satisfied) 

For every coefficient plane, set list of significant pixels (LSPipiane) as NULL, then add the 
coordinate (i,j) € H plane to LIPipiane and those that have descendants to LISipiane, at the same 
time, set their type to be A. 

2. Sorting pass 

If n. max; <n, only sort C4(i, j) just as conventional SPIHT algorithm. 

If n. max; «n«n, maxs, sort C4(i,j) and then sort C;(i,j) in succession. 

If n «n, max, , sort C,(i,j) and then sort Cy(i, j) , C4(;, j) in succession. 

3. Refinement pass 

If n. max; <n, for any entry (ij) in LSP, except those included in the last sorting pass, 
output the nth most significant bit of C,(i, j) . 

If n. max; «n «n, maxs, for any entry (ij) in LSP, and LSP; except those included in the 
last sorting pass, output the nth most significant bit of C,(i,7) and C}(i, j). 

If n «n. max; , for any entry (ij) in LSP,,LSP, and LSP, except those included in the last 
sorting pass, output the nth most significant bit of C,(i,7), C,(i,j) and C4(i,j) . 

4. Quantization updating 

Decrement n by 1 and go to step 2 until acquire the desired compression ratio. 

The code stream of embedded 3D-SPIHT encoding can be illustrated as figure.2. The 
decoding process and encoding process are just the same. So, it's easy to present the 
decoding process according to the 3D-SPIHT encoding process. The code stream of 3D- 
SPIHT is more compact and efficient for the three mixed coefficient planes are interleaved 
during the encoding process. 


n=n-| 


C; C5 C3C; C5 C3... 





n-n max, n max, n max, 1 


Fig. 2. Code stream illustration for embedded 3D-SPIHT. 
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Additionally, it is necessary to be mentioned that there have been two kinds of 3D-SPIHT 
algorithms before. One is proposed for video compression by extending the conventional 
SPIHT algorithm to 3D case directly and encoding the 3D-DWT wavelet coefficients of video 
data, so the SOF is defined as 8 splitting [7]. The other is proposed for compression of 
multispectral images, which make some amendments of the conventional SPIHT by adding 
one spectral child to every baseband coefficient, so its SOF is still 4 splitting [8]. The 
proposed 3D-SPIHT embedded coding in this book is very different from the two existing 
3D-SPIHT algorithms, which encodes 1 or 2 or 3 coefficient planes of the 3 mixed coefficient 
plane sequentially by adopting 3 independent thresholds. 


4. Unsymmetrical SPIHT codec 


It is easy to find that the set partitioning is keeping on 4 partitioning for conventional SPIHT 
codec, that is to say the same set partitioning in vertical and horizontal direction. This means 
that the wavelet decomposition only can be implemented under the constraint of the same 
decomposition level in vertical and horizontal direction, which will affect the redundancy 
removing of image data efficiently and the final compression performance. The 
unsymmetrical SPIHT codec, adopting set 2-partitioning or set 4-partitioning according the 
requirements, doesn’t require the same decomposition level in vertical and horizontal 
direction. So, DWT can be implemented with each highest feasible decomposition level in 
vertical and horizontal direction respectively and then the spatial redundancy can be 
removed efficiently. 

Because the flow chart of unsymmetrical SPIHT codec is very near to conventional SPIHT 
codec, only the difference is given. 

Assume the image dimensions are W*H and H=h* 2HLed W = w*2Weel here 
HLevel,WLevel are the highest feasible decomposition level in vertical and horizontal 
direction respectively. Usually, there exists HLevel#WLevel when the image size are 
unsymmetrical. For the conventional SPIHT, the highest feasible decomposition level only 
can be the lower one of HLeveland WLevel, so, Level = min(WLevel, HLevel). But for 
unsymmetrical SPIHT codec, the highest wavelet decomposition levels in vertical and 
horizontal direction are HLevel and WLevel respectively. 

All set partitions for conventional SPIHT are set 4 partitioning just as the formula 2, but the 
set partitions for unsymmetrical SPIHT are set 2-partitioning or set 4-partitioning, just as the 
following formula: 


ci2j2;-1) ifHy <i <2 


O(i,j) = c(2i,jc(2i-1,j) if iWysj« = (18) 
c(2i,2j),c(2i+1,2j),c(2i,2j+1),c(2i+1,2j+1) otherwise 








H W 


Here, Ho = pitenel-Wieoel+T ’ 0 = o WLevel-HLevel 1 ` 


When HLevel = WLevel , the unsymmetrical SPIHT codec will completely degenerate into 
conventional SPIHT codec. When HLevel -1 or WLevel=1, only set 2-partitioning are 
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implemented in one direction, the coefficients can be encoded line by line independently. 
The wavelet transform in another direction is only be used as compacting the image energy. 
Fig 3 and Fig 4 are the illustrations of conventional SPIHT encoding and unsymmetrical 
SPIHT encoding at HLevel=3 or WLevel=5. 





Fig. 3. The illustration of set partitioning using conventional SPIHT encoding for 
unsymmetrical image size 


MIIIIIIIE————————ÉMNIIIIITIITITITIT] 


MIIIIIIIIIIIIIIII 





Fig. 4. The illustration of set partitioning with unsymmetrical SPIHT encoding for 
unsymmetrical image size 
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5. 1D SPIHT codec 


In section 4, unsymmetrical SPIHT codec is detailed described, which also need 2D image 
data or image block. But, in real time image transmission or scan display system, the image 
data are usually transmitted or displayed line by line. In order to use conventional SPIHT or 
unsymmetrical SPIHT, it needs many line buffers to store the previous image data. In 
hardware, it is a high burden for the costly RAM. So, 1 line image data compression method 
will have the precedence over other block based compression methods, such as the 1D DWT 
followed 1D SPIHT codec which will be addressed in the following. 

After 1D DWT, the wavelet coefficient also has the natural pyramid characteristic: every 
pixel of the high frequency subband has its 2 corresponding pixels in its adjacent level high 
frequency subbands in position, which means that only set 2-partitioning can adopted. The 
illustration is given in fig3. 


| | ee 


Fig. 5. The illustration of set partitioning with 1D SPIHT encoding for 1 line image. 


The SOT and set partitioning can be written as formula 19 and 20. 


D(i) = O(i) + L(i) (19) 
Lü)-YD(k) (k)eO(i) 
O(i) = {c(2i), c (2i+1)} (20) 


From fig 5 and formula 19 and 20, we can see that 1D SPIHT use set 2 partitioning to encode 
1D DWT coefficients, which only need 1 line buffer RAM but leave another dimensional 
redundancy un-removed. 


6. Conclusion 


In this chapter, SPIHT and it’s derivatives or its modification methods, such as 3D-SPIHT, 
3D-SPIHT Embedded method, Unsymmetrical SPIHT and 1D-SPIHT are described, 
which can overcome the disadvantages of traditional SPIHT codec and meet the 
specific requirements for the real applications. Fig.6 gives the derivative relationship of 
traditional SPIHT and its several modified methods. We can see that SPIHT is the 
foundation for traditional symmetrical image, but it can’t meet the requirements at some 
specific applications, such as multi-channel image, strip image (unsymmetrical image), 
even image line. But, its modified version can meet some specific requirements in real 
applications. 
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Methods Dimensions Application case 
3D-SPIHT 3D Video 
3D-SPIHT 2.5D' Multi-Channel images 
Embedded 
SPIHT 2D Image 
Unsymmetrical 1.5D' Unsymmetrical image 
SPIHT /Strip image 
1D-SPIHT 1D 1 Line image 


Or other 1D signal 


* mean the corresponding dimension is limited compared other dimensions 


Fig. 6. Family of SPIHT and its derivatives. 


7. References 


[1] Said A. and Pearlman W.A., A new fast and efficient image codec based on set 
partitioning in hierarchal trees, IEEE Trans. on Circuits and Systems for Video 
Technology, 1996, 6(3): 243-250. 

[2] D. E. Dudgen, R. M. Mersereau, Multidimensional digital signal processing, New 
Jersey : prentice-Hall, Inc.,1984. 

[3] Zhu Yangiu > Chen Hexin » Dai Yisong. Compression Coding of color image via 3-D 
matrix transform [J] ACTA EI ECTRONICA SINICA > 1997 > 25(7):16-21. 

[4] Y. Shoham, A. Gersho, Efficient bit allocation for an arbitrary set of quantizers, IEEE 
Transaction on Acoustics,Speech, and Signal Processing, 1988,36(9):1445-1453 

[5] P. Prandoni, M. Vetterli, R/D optimal linear predication, IEEE Transaction on Speech 
and Audio Processing, 2000,8(6):646-655. 

[6] A. Aminlou, O. Fatemi, Very fast bit allocation algorithm, based on simplified R-D curve 
modeling, IEEE ICECS 2003, 112-115. 

[7] B. J. Kim, Z. Xiong, W. A. Pearlman,Low bit-rate scalable video coding with 3-D set 
partitioning in hierarchical trees(3-D SPIHT), IEEE Transaction on Circuilt and 
System for Video technology, 2000, 10(8): 1374-1387. 

[8] P. L. Dragotti, G. Poggi and A. R. P. Ragozini, Compression of multispectral images by 
three-dimensional SPIHT algorithm, IEEE Transaction on Geoscience and Remote 
sensing, 2000, 38(1):416-428. 


www.intechopen.com 


78 Discrete Wavelet Transforms: Algorithms and Applications 





[9] W.-C. Zhang, Y.-F. Wang, G.-H. Hu,.Compression of multi-polarimetric SAR intensity 
images based on 3D-matrix transform, IET- Image processing, 2008,2(4):194:202. 

[10] Zhang Zhi-hui, Zhang Jun, Unsymmetrical SPIHT Codec and 1D SPIHT Codec, 
International Conference on Electrical and Control Engineering (ICECE), 2010, 
wuhan 2498:2501. 


www.intechopen.com 


Discrete Wavelet Transforms - Algorithms and Applications 
Edited by Prof. Hannu Olkkonen 


ISBN 978-953-307-482-5 

Hard cover, 296 pages 

Publisher InTech 

Published online 29, August, 2011 
Published in print edition August, 2011 





The discrete wavelet transform (DWT) algorithms have a firm position in processing of signals in several areas 
of research and industry. As DWT provides both octave-scale frequency and spatial timing of the analyzed 
signal, it is constantly used to solve and treat more and more advanced problems. The present book: Discrete 
Wavelet Transforms: Algorithms and Applications reviews the recent progress in discrete wavelet transform 
algorithms and applications. The book covers a wide range of methods (e.g. lifting, shift invariance, multi-scale 
analysis) for constructing DWTs. The book chapters are organized into four major parts. Part | describes the 
progress in hardware implementations of the DWT algorithms. Applications include multitone modulation for 
ADSL and equalization techniques, a scalable architecture for FPGA-implementation, lifting based algorithm 
for VLSI implementation, comparison between DWT and FFT based OFDM and modified SPIHT codec. Part II 
addresses image processing algorithms such as multiresolution approach for edge detection, low bit rate 
image compression, low complexity implementation of CQF wavelets and compression of multi-component 
images. Part III focuses watermaking DWT algorithms. Finally, Part IV describes shift invariant DWTs, DC 
lossless property, DWT based analysis and estimation of colored noise and an application of the wavelet 
Galerkin method. The chapters of the present book consist of both tutorial and highly advanced material. 
Therefore, the book is intended to be a reference text for graduate students and researchers to obtain state- 
of-the-art knowledge on specific applications. 


How to reference 
In order to correctly reference this scholarly work, feel free to copy and paste the following: 


Wenchao Zhang (2011). Several Kinds of Modified SPIHT Codec, Discrete Wavelet Transforms - Algorithms 
and Applications, Prof. Hannu Olkkonen (Ed.), ISBN: 978-953-307-482-5, InTech, Available from: 
http://www. intechopen.com/books/discrete-wavelet-transforms-algorithms-and-applications/several-kinds-of- 
modified-spiht-codec 


INTECH 


open science | open minds 


InTech Europe InTech China 

University Campus STeP Ri Unit 405, Office Block, Hotel Equatorial Shanghai 

Slavka Krautzeka 83/A No.65, Yan An Road (West), Shanghai, 200040, China 

51000 Rijeka, Croatia FESTER R65 LEER RER AT EN A405 TE 
Phone: +385 (51) 770 447 Phone: +86-21-62489820 


www.intechopen.com 


Fax: +385 (51) 686 166 Fax: +86-21-62489821 
www.intechopen.com 


© 2011 The Author(s). Licensee IntechOpen. This chapter is distributed 
under the terms of the Creative Commons Attribution-NonCommercial- 
ShareAlike-3.0 License, which permits use, distribution and reproduction for 
non-commercial purposes, provided the original is properly cited and 
derivative works building on this content are distributed under the same 
license. 


